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Abstract 

We derive high-order terms in the asymptotic expansions of the steady- 
state voltage potentials in the presence of a finite number of diametrically small 
inhomogeneities with conductivities different from the background conductiv¬ 
ity. Our derivation is rigorous, and based on layer potential techniques. The 
asymptotic expansions in this paper are valid for inhomogeneities with Lipschitz 
boundaries and those with extreme conductivities. 
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1 Introduction 

Let U be a bounded domain in M'^, d > 2, with a connected Lipschitz boundary 
dU. Let v denote the unit outward normal to dU. Suppose that U contains a finite 
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number m of small inhomogeneities each of the form Di = zi + eBi, where 

Bi, I = 1,... ,m, is a bounded Lipschitz domain in containing the origin. We 
assume that the domains are separated apart from each other and apart 

from the boundary. More precisely, we assume that there exists a constant cq > 0 
such that 


(1.1) \zi — zi>\ > 2co > 0,V / / /' and dist(2;z, 911) > 2co > 0,V/, 


and e, the common order of magnitude of the diameters of the inhomogeneities, is 
sufficiently small, that these inhomogeneities are disjoint, and that their distance to 
W^\Q is larger than cq. We also assume that the ’’background” is homogeneous with 
conductivity 1 and the inhomogeneities Di have conductivities ki, ki ^ 1, 1 < I < m. 

Let Ue denote the steady-state voltage potential in the presence of the conduc¬ 
tivity inhomogeneities, i.e., the solution to 


( 1 . 2 ) 


/ mm 

V- (x(ll \ U ^) + E kixiDi)) Vu, = 0 

< i=i 1=1 


due I 

> du ^ 


in n. 


Let U denote the ’’background” potential, that is, the solution to 


(1.3) 


AC/ = 0 in n, 
, du ^ 


The function g represents the applied boundary current; it belongs to Lq((911) = {5 € 
L^(dQ), Jqq g = 0}. The potentials, Ue and U, are normalized by Ue = C/ = 0. 

The main achievement of this paper is a rigorous derivation, based on layer 
potential techniques, of high-order terms in the asymptotic expansion of Uelao as 
e ^ 0. The leading order term in this asymptotic formula has been derived by 
Cedio-Fengya et al. ; see also the prior work of Friedman and Vogelius |]lj] for the 
case of perfectly conducting or insulating inhomogeneities. The main result of this 
paper is the following full asymptotic expansion of the solution for the case m = 1 : 


Theorem 1.1 Suppose that the inhomogeneity consist of single component and let 
Us be the solution of The following pointwise asymptotic expansion on dQ 
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holds for d = 2,3: 


n n—\i\+l 




(1.4) 


Ue{x) 

14=1 lil=i 

n-\-2—\i\ — \j\—d 


1+ Y. (^^^-^Qp){d^U{z))]Mi,diN{x,z) 


p=i 


+0{e^+^), 


where the remainder is dominated by Ce‘^^^\\g\\i 2 (^g^'^ for some C indepen¬ 

dent of X ^ dfl. Here N{x,z) is the Neumann function, Mij, i,j G are the 
generalized polarization tensors defined in and the operator Qp is defined in 

(EH)- 

We have a similar expansion for the solutions of the Dirichlet problem (Theorem 


The derivation of the asymptotic expansions for any fixed number m of well 
separated inhomogeneities (these are a hxed distance apart) follows by iteration of 
the arguments that we will present for the case m = 1. In other words, we may 
develop asymptotic formulas involving the difference between the fields and U on 
dit with I inhomogeneities and those with I — 1 inhomogeneities, I = m,... , 1 , and 
then at the end essentially form the sum of these m formulas (the reference fields 
change, but that may easily remedied). The derivation of each of the m formulas is 
virtually identical. 

We also note that the asymptotic expansion o is valid for inhomogeneities 
with zero or inhnity conductivity (cavity or perfect conductor). Precise definitions 
of generalized polarization tensors (GPT) associated with the domains Bi and the 
conductivities ki will be given at the end of section 3. These GPT seem to be natural 
generalizations of the tensors that have been introduced by Schiffer and Szego ||2^ 
and thoroughly studied by many other authors 11 , H], 0 , i- (See section 3.) 

The higher-order terms are essential when VU{zi) = 0, for then the leading 
order term in the asymptotic expansion of rte|an, given in ||^], vanishes. We remind 
the reader that, for general current inputs 5 , Vf7 vanishes at some ’’critical points” 
inside 11 . 

The proof of our asymptotic expansion is radically different from the ones in 
a, §, and |@. It is based on layer potential techniques and a decomposition 
formula of the steady-state voltage potential into a harmonic part and a refraction 
part. This formula is due to Kang and Seo |T^. What makes our proof particu¬ 
larly original and elegant is that the rigorous derivation of high-order terms follows 
almost immediately. The extension of the techniques used in [1T|, 0 , and @1 to 
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construct higher-order terms in the expansion of Uf\dQ, as e —> 0 seems to be labo¬ 
rious. Furthermore, the general approach developed in this paper could be carried 
out to obtain more precise asymptotic formulas for the full Maxwell’s equations and 
for the equations of linear elasticity than those derived in Q and |j|]. The method 
of this paper also enables us to extend the asymptotic expansions to the cases of 
inhomogeneities with Lipschitz boundaries. Previously, the leading order term was 
derived under the assumption that inhomogeneities are (7^’“ smooth |l^ |Q]. We 
note that our method works as well even when the inhomogeneities have extreme 
conductivities (A: = 0 or A: = oo). 

Let us now explain what makes this asymptotic formula interesting in the elec¬ 
trical impedance tomography (FIT). It is well-known that the ultimate objective of 
FIT is to recover, most efficiently and accurately, the conductivity distribution in¬ 
side a body from measurements of current flows and voltages on the body’s surface. 
The vast and growing literature reflects the many possible applications of FIT, e.g. 
for medical diagnosis or nondestructive evaluation of materials |6|. In its most gen¬ 
eral form FIT is severely ill-posed and nonlinear. Taking advantage of the smallness 
of the inhomogeneities, Cedio-Fengya et al. |@] have used the leading order term 
in the asymptotic expansion of u^\qq, to find the locations zi, I = 1 ,. .. ,m of the 
inhomogeneities and certain properties of the domains Bi, I = 1 ,... ,m (relative 
size, orientation). The algorithm proposed by Cedio-Fengya et al. 0 is based on 
a least-square algorithm. Ammari et al. have also utilized this leading order 
term to design a variationally based direct reconstruction method. The new idea in 
1 ^ is to form the integral of the ’’measured boundary data” against harmonic test 
functions and choose the input current g so as to obtain expression involving the in¬ 
verse Fourier transform of distributions supported at the locations zi, I = 1,... ,m. 
Applying a direct Fourier transform to this data then pins down the locations. This 
approach is similar to the ideas used by Calderon in his proof of uniqueness of 
the linearized conductivity problem and later, by Sylvester and Uhlmann in their 
important work |^] on uniqueness of the three-dimensional inverse conductivity 
problem. An other algorithm that makes use of an asymptotic expansion of the 
voltage potentials has been derived by Briihl et al. Q. This algorithm is in the 
spirit of the linear sampling method of Colton and Kirsch |^. 

In all of these algorithms, the locations zi, I = 1,... ,m of the inhomogeneities 
are found with an error 0(e) and little about the domains Bi can be reconstructed. 
Making use of higher-order terms in the asymptotic expansion of u^ldn we certainly 
would be able to reconstruct the small inhomogeneities with higher resolution from 
boundary information about specific solutions to O- Perhaps, more importantly, 
this would allow us to identify quite general conductivity inhomogeneities without 
restrictions on their sizes. 

The use of higher-order terms in the asymptotic expansion of Ue\dn may also be 
decisive in dramatically improving the algorithm of Kwon et al. , that is based on 
the observation of the pattern of a simple weighted combination of an input current 
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g of the form g = a ■ v for some constant vector a and the corresponding output 
voltage. We also believe that the use of such higher-order terms would improve the 
algorithm of Mast et ah that uses eigenfunctions of the scattering operator. 

This paper is organized as follows. In section 2, we collect some notations and 
preliminary regarding layer potentials. In section 3, we introduce the generalized 
polarization tensors associated with the domains Di and the conductivities ki. In 
section 4, we provide a rigorous derivation of high-order terms in the asymptotic 
expansion of the output voltage potentials. For reasons of brevity we restrict a 
significant part of this derivation to the case of a single inhomogeneity (m = 1). 
The proof in the case of multiple well-separated inhomogeneities may be derived by 
a fairly straightforward iteration of the arguments we present, however we leave the 
details to the reader. 


2 Layer Potentials for the Laplacian 


Let us first review some well-known properties of the layer potentials for the Lapla¬ 
cian and prove some useful identities. 

The theory of layer potentials has been developed in relation to the boundary 
value problems. Let D be a bounded domain in d > 2. We assume that dD is 
Lipschitz. Let r(x) be the fundamental solution of the Laplacian A: 


( 2 . 1 ) 


r(x) = 


^InN, 


, (2 - d)oJd 


\2-d 


d = 2, 
d > 3, 


where ujd is the area of (d — 1) dimensional unit sphere. The single and double layer 
potentials of the density function on ZD is defined by 


(2.2) 

Sd4>{x) 

:= [ r{x - y)4>{y)da{y), xGM^, 



JdD 

(2.3) 

Dnfpix) 

r o 

:= / —T{x-y)(t){y)da{y), x&W^\dD 

JdD dvy 


For a function u defined on R'^ \ dD, we denote 


, u{x) := lim (Vu(x ± tUx)-, ^x)-, x € dD 
du^ t^o+ 

if the limit exists. Here is the outward unit normal to dil at x. 
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The proof of the following trace formula can be found in |11, |^, (for Lipschitz 
domains, see |^|): 


(2.4) ^^DcPix) = i±h + IChmx), 

(2.5) (^<(>)l± = {T^I + ICd)4>{x), X £ dD, 

where 

ICofpix) = —p.v. f —^f^(/>(y)dcr(y) 

JdD \x - yr 

and is the L^-adjoint of JCd- When dD is Lipschitz, ICd is a singular integral 
operator and known to be bounded on L‘^{dQ) p. Let LQ{dD) := {/ £ L‘^{dD) : 
JdD ~ following results are due to Verchota and Escauriaza-Fabes- 

Verchota. 


Theorem 2.1 ( [lTo| , |]2^ ) XI — K.*j^ is invertible on L^{dD) if |A| > and for 
A £ (—oo, —i] U (i, oo), XI — JCf) is invertible on L'^{dD). 


For proofs when dD is smooth, see [pT[|, ||I3|. 


The following theorem was proved in 


0 , 17 ]. 


Theorem 2.2 Suppose that D is a domain eompaetly contained in Q with a con¬ 
nected Lipschitz boundary and conductivity k. Then the solution u of the problem 


( V • ((1 + (A: - l)x(D))Vu) =0 in Q, 

(^•6) I du, 

V di' ^ 

is represented as 

(2.7) u{x) = II{x)-\-SD(j>{x)^ X £ n, 
where the harmonic function H is given by 

(2.8) H{x) =-Su{g){x)+Vu{f){x), x £ fl, /:=u|an, 
and 4> £ L‘f^{dD) satisfies the integral equation 

k 1 dH 

£3) <2(rvij^-':BW = ^U »»««• 

Moreover, V n £ N, there exists a constant Cn = C(n,fl, dist(D,dQ)) independent 
of \D\ and k such that 


( 2 . 10 ) 


.^llc"(D) — ^n\\g\\L^(dn)- 


6 








Proof. The representation formula (^) was proved in |]^, |^. ( p.lO|) was 

proved in |^] for d = 2 and it is easily seen that the same proof works for d = 3. 
We only need to check carefully whether the constant Cn in the estimate (|2.10| ) is 
independent of \D\. Before doing this, let us point out that the harmonic function H 
can be computed explicitly from the boundary measurements {^\q^i'u\qq) and the 
density (p is uniquely and explicitly determined by the domain D and the harmonic 
function H. The decomposition of the function u into a harmonic part H and a 
refraction part Sjjfp unique |0, 171. The representation formula (2.7) seems to 


inherit geometric properties of D. 

Suppose that dist(Zl, 90)) > 2co for some constant cq > 0. From the definition 


of H in (|2.8|) it is easy to see that 


( 2 . 11 ) 


|77||^n(£)) < Cn 


L^(aQ) + Ikiso 


where depends only on n, 90, and cq. Let a be a vector field supported in the 
set dist(x, 90) < 2co such that a ■ u{x) > 5 for some (5 > 0 for all x G 90. Using the 
Rellich identity with this a, we can show that 

\\-^\\L^{dn) < C'(ll5llL2(0n) + 

where C depends only on 90 and cq and T{x) is the tangent vector to 90 at x. (See 


the proof of Lemma 2.1 of [12| for the detail of the proof.) Observe that 

^ ^ + ik- l)x(Zl))Vu • Vu dx 

I qu da 

Jan 




= C 


< C\\g\\L2{an)\W\an\\L'2(an)- 

Since u da = 0, if follows from the Poincare inequality that 


du 


I 190 111,2 (QQ) < Ull —11^2(9^) 


Thus we obtain 




and hence 


^l9o||L2(an) < C'llfi'llL2(an)- 


From (2.11) we 


finally obtain ( p.lO|) . 


□ 
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Using above representation we can derive a formula similar to (1.4) which is 
potentially useful in detecting the inhomogeneities (see the remark at the end of 
this paper). However it uses the function H which depends on D and hence on 


e. Thus in order to derive (^) we will transform it to representations using only 
background potentials. 

Let N{x, y) be the Neumann function for A in H corresponding to a Dirac mass 
at z. That is, N is the solution to 

A.xN{x,z) = —5z in H, 
dN. _ 1 

In addition, we assume that 

(2.12) / N{x,y)da{x) = 0 for y G H. 

Jon 

Let us hx one more notation: For D, a subset of let 


NDf{x):= N{x,y)f{y)da{y). 

Jan 

The following lemma relates the fundamental solution with the Neumann func¬ 
tion. 

Lemma 2.3 For z G D and x G d^l, let r 2 (x) := r(x — z) and Nz{x) := N{x,z). 
Then 

(2.13) + ^n){d^z)ix) = r 2 (x) modulo constants, x G dfl, 

or to be more preeise, for any simply connected Lipsehitz domain D eompactly eon- 
tained in D and for any g G L‘^{dD), we have 

(2.14) f {-l-I + }Cn){Nz){x)g{z)da{z) = f Tz{x)g{z)da{z), Vx G 

JdD ^ JdD 


Proof. Let / G ^^(SH) and define 


u{z) := {{--I + ICn){Nz),f)an zeTl. 


u{z)=[ N{x,z){-^I+ K:^)f{z)da{z). 
Jan z 


Then 






Therefore, Au = 0 in and ^h)f- Thus by (2.4) we have 


u{z) — Snf{z) = constant, z € Cl. 


Thus if 5 G LQ{dCl), then we obtain 


/ / + l<^n){Nz)ix)g{z)da{z)f{x)da{x) 

Jan JdD ^ 

= [ [ C^{x)g{z)da{z)f{x)da{x). 

Jan JaD 


Since / is arbitrary, we have (2.13) or equivalently, ( |2.14| ). This completes the proof. 

□ 

Let g € Lg(5n). Let U{y) := Jq^ N{x,y)g{x)da{x). Then U satisfies 

Af7 = 0 in n. 


(2.15) 


~^\dn = g ^ To(9fl), 

/ U {x)da{x) = 0. 

Jan 


Theorem 2.4 The solution u of ( \2. dj ) can be represented as 
(2.16) u{x) = U{x) — Nd(P{x), X € dCl 

where f is defined in (|A ^). 


Proof. By substituting (|2.7D into the equation (|2.8|) , we obtain 

H{x) = -Sn{g){x) + VniHlon + (5z)(?i)|af2)(a:), x e Cl. 
It then follows from (|2.5|) that 


(2-17) {-I - K-n){H\QQ,) =-{Sng)\9n + {-I + }Cn){{SD4>)\an), on 911. 

Since U = -Sn{g) + PniUlan) in Cl, we have 


(2.18) 


- }(^n){U\an) = -{Sng)\9n. 


Since (p G Lq^OD), it follows from (2.12) that 


(2.19) 


~ ^^){{^D4>)\an) = {SDfi)\an- 
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From ( 2.17|) , (p.l 8 ), and ( 2.19 ), we conclude that 


- ^Q){H\dQ, - U\qq, + (-/ + D4>)\dQ)) = 0. 


Therefore, we have 


(2.20) H\qq, - U\dQ, + (-/ + 1Cq)[{N£,( t))\QQ) = C (constant). 

Note that {\l + K,n){{NDct>)\dn) = (^^D</>)|an + {SD4>)\dn- Thus we get from ( ^ 
and (|] 2 ^ 


(2.21) u\gn = U\qq - {ND4>)\dn + C. 


Since all the functions entering in 
and the theorem is proved. 


(| 2 . 21 ) belong to ^^(90), we conclude that (7 = 0, 

□ 


We have a similar representation for solutions of the Dirichlet problem. Let 
G{x,y) be the Green function for the Dirichlet problem, i.e., the function V defined 
by V{x) := af fa) y)f{y)do'{y) is the solution of the problem AV = 0 in D and 

= / for any / G L^(5D). Then we have the following representation theorem. 


Theorem 2.5 


( 2 . 22 ) 




dr,{y) 

dv{y) 


){x) 


dG, 

dv{x) 


(a:), 


X G 9D, z G D. 


Let u be the solution of ( [^.6[ ) with the Neumann condition replaced by the Dirichlet 
condition ul^n = /. Then u can he represented as 


(2.23) 



dV 

dv 


{x) - GdH^), 


X G dLl 


where (j) is defined in 


and 


GD(t){x) := [ 

JdD 


dG 

dn{y) 


{x,y)4>{y)dcr{y). 


Theorem |2.5| can be proved in the same way as Theorem In fact, it is simpler 
because of the solvability of the Dirichlet problem, or equivalently, the invertibility 
of (^/ + JCq). So we omit the proof. 


3 Generalized Polarization Tensors 

In this section we introduce the generalized polarization tensors (GPT) associated 
with a domain B and a conductivity k. These GPT are the basic building block for 
the asymptotic expansions in this paper. 
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Let -B be a Lipschitz bounded domain in and the conductivity of B be fc 
{k 7 ^ 1). The polarization tensor is M = {rriij), I < i, j < d, is dehned by 

f d 

^ij\B\+ yi^—il)j{y)da{y) , 

JdB 

where is the unique solution of the following transmission problem: 

' A'i/’j(x) = 0, xGBUM'^\B, 

'4^j\+ — V'jl- = 0 on dB, 

= i^j ondB, 

—> 0 as \x\ —> oo. 


rriij ■= {k - 1) 


See [|^], 0, and |Q]. One can easily check using 


that 




Using (^) again, we have 


= - / yi(b - h’bKh - ic'BrHyj)(y)'rr{y) 

= -f yii^jdaiy) + {X - [ yi{XI - }C*B)~^{iyj){y)da{y) 

Job ^ JdB 

= -6ij\B\ + f yi{XI - JC%)~\uj){y)da{y). 

K — i- JdB 

Therefore we prove that the polarization tensor M associated with B and k is given 


(3.1) rriij = [ yiiXI-KTs) ^{vj){y)da{y). 

JdB 

Recall A :=^. 

For a multi-index i = (r, • • • ,id) ^ let 9*/ = f and x* := ■ 

For i,j G N'^, we dehne the generalized polarization tensor Mij by 

(3.2) Mij := [ y^4>i{y)da{y), 

JdB 

where (pi is defined by 

pi{x) := {XI - lC%)~^{^Uy ■ Vy*)(x), x G dB. 
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4 Derivation of the Full Asymptotic Formula 


In this section we derive our asymptotic formula (p..4|). As stated in the introduction, 
we restrict our derivation to the case of a single inhomogeneity {m = 1). We only give 
the details when considering the difference between the fields corresponding to one 
and zero inhomogeneities. In order to further simplify notation we assume that the 
single inhomogeneity D has the form D = eB + z, where z € and i? is a bounded 
Lipschitz domain in containing the origin. Suppose that the conductivity of D 
is k. Let A := 2 {k-i) ■ Then by (^) and (|2.9D , the solution u of (|2.6|) takes the form 

BH 

u{x) = U{x) - Nd{\I - |aD)(a:), x G 89., 

where U is the background potential given in (|2.15 ). 

Dehne 

1 


Hn{x) ■.= Y,-{d^H){z){x-z)\ 


j =0 


Then we have from ( 2.IC ) that 

dH dH, 


du dv 


Note that 
(4.1) 


if 


< sup \VH{x)-VHn{x)\\dD\^l‘^ 
L'2{dD) x&dD 

<C\\gU^an)^^\dD\^/\ 

f hda = 0, then f {\I — lCy)~^hda = Q. 

IdD JdD 


If fgjy hda = 0, then we have for x G 89 that 

\Nd{XI - ICh)-^h{x)\ = [ [N{x-y)-N{x-z)]{XI-JCh)-^hiy)daiy) 

JdD 

<Ce\8D\^/^\\h\\L,^Q^y 


It then follows that 


B M 

Nd{XI - {-^\oD 


Therefore, we have 


dH^ 

8u 


\dD 


)ix] 


< Ce\8D\^/^ 

< C\\g\\L^(dn)<^' 


8H 8H, 
8u 

.d-\-n 


8v 


L^{dD) 


B M 

(4.2) u(x) = U{x) - Nd{XI - )Ch)-H^\dD){x) + 0(6'"+") 


X ^ 89 
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where 0(e'^+”) term is dominated by C'|| 5 ||L 2 (ar 2 )e‘^^"' for some C depending only on 
cq. Note that 

rtf-f ^ 1 

(A/ - K*j,)-\^\aD){x) = Y,{d^H){z){XI - ICh)-\-^ ■ V(x - zy)ix). 

|i|=i 


Since D = eB + z, one can prove by using the change of variables y = that 
(A/ - • V(x - zyyx) = el^i-^A/ - lCj,Y\\vy ■ - z)). 

Put 

(4.3) := (A/-■ W)(a;), x^dB. 

Then we get 


(4.4) 


ND{\I-K.l)-\^\aD){x) 


^(a*i/)(z)el*l M N{x,y)(j)i{e ^{y-z))da{y) 


|i|=i 

n 


IdD 


^(a*i/)( 2 ;)el*l+'^ ^ [ N{x,ey + z)(t>i{y)da{y). 
H=i 


We now expand N{x,ey + z) asymptotically as e ^ 0. By (|2.13D we have the 
following relation: 

{—-I + K-Yj[X^{■ ^ ^y + z)]{x) = T{x — z — ey) modulo constants, x ^ dVL. 


Since 


H-OO 


1 


r(x-ey) = ^ -e>-’la^(r(x))y^ 


we obtain 



zW 


= {-U + lCn) 



(x). 


Since N{x, w)da{x) = 0 for all tc G fl, we have the following asymptotic expan¬ 
sion of the Neumann function which is of independent interest. 


13 







Lemma 4.1 For x G dil., z & il., and y G dB, and e —> 0, 


(4.5) 


+ 00 ^ 

N{x, ey + z) = E —e^^^diN{x, z)yF 
\j\=o^' 


We now have from (4.4) 


f) M 

NoiXI - ICh)-\^\9D){x) 


+ CXD 


1 


^(9*i/)(2;)el*l+'^ 2 -e^^^dlN{x,z) / y^(j)i{y)da{y). 

1 * 1=1 lj|=o-^' 


Observe that 


2 i(8‘i?)(2)A(!,*) = A J 


= 0 , 


HH ^ \l*l=; 

and therefore, by the Green’s theorem, it follows that 

1 


Yz:{d^H){z)V{y^)-u{y)da{y) = 0. 

\i\=i *• 

Thus, in view of (^), the following identity holds by using observation o 

(4.6) V(5*iJ)(z) / 4>i{y)da{y) = 0, V 

nr, JOB 


i > 1 . 


\i\=l 


In fact, this follows immediately from (4.1). Recall now that Mjj = / y^(j)i{y)da{y) 

JdB 

is the generalized polarization tensor associated with the domain B and the conduc¬ 
tivity k to obtain the following pointwise asymptotic formula: for x G dVt, 


72—121 + 1 


(4.7) u{x) = U{x)-e'^-‘^Y \(^"^^^^\d"H){z)MijCp,N{x,z)+0{e'^+'^). 

14=1 |ii=i 


Observing that the formula (4.7) still contains d^H factors, the remaining task 
is to convert (4.7) to a formula given solely by U and its derivatives. 

As a simpliest case, let us now take n = 1 to find the leading order term in the 
asymptotic expansion of u|gn as e —> 0. From G) and (|2.16|) , we get 


+ - ^+“(90) < 11'/’ + 2 ( 925 ) < C'e'^||fl'llL2(9n) 
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for some C depending only on Q and cq. It then follows from the maximum principle 
that 

\\H - C/||loo(q) < Ce^\\g\\L2(^QQy 

Then, from the mean value property of harmonic functions, we obtain 

\VH{z)-VUiz)\<Ce%h2^9^y 

It thus follows from ( [4.7D that 

(4.8) u{x) = U{x) {d^U){z)Mijd^N{x, z) + 0(e‘^+^), x G dn, 

N=blil=i 

which is in view of exactly the formula derived in [14| and when D has (7^’“ 
boundary. 


We now return to (^). Substitution of (O) into (|2.8|) yields that, for any 
X G n. 


(4.9) 


n—1«| + 1 


1 u. 


H{x) = U{x) - e'^-2 Y -e^^^^^^Kd"H){z)MijVn{diN{-, z)){x) + 0(e‘^+"). 
14=1 lil=i 

In (4.9) the remainder 0(e'^+”') is uniform in the norm on any compact subset of 
n for any n and therefore 

n n-|il+l 

(4.10) {d^H){z) + Y^''~^ Y. e^'^^^^Kd'H){z)P^y = {d^U){z)+0{e^+^), 

14=1 lil=i 


for all I G with \l\ < n where 

(4.11) P,ji = ^Mi,dlVn{diN{; z))U=, 

Dehne the operator 


■ (ii/)zeN‘*,|4<n + e 


d-2 


n n—|i| + l 

E E ' 

14=1 lil=i 


14+lil 


ViP 


ijl 


leN‘i,\l\<n 


Observe that 

V, = I + e^Vi + ... + 

Defining Qp,p = 1,... ,re - 1, by 

(4.12) 

(/ + e^Vi + ... + e”+''-^iP„_i)-^ = 1 + e^Qi + ... + e^+^-^Qn-i + 
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we finally obtain that 


(4.13) ((a'if){4))ia.,|,|<„ = (/ + E E''+’’-‘Cp)({8‘t')(4))i6H-,|.|<„ + O)!"”) 


P=l 


which yields the main result of this paper stated in Theorem 1.1. 

We also have a complete asymptotic expansion of the solutions of the Dirichlet 
problem: 


Theorem 4.2 Suppose that the inhomogeneity consist of single component and let 
u be the solution of with the Neumann condition replaced by the Dirichlet 

condition u|an = /. Let V be the solution of AV = 0 in Q with V\qq, = /. The 
following pointwise asymptotic expansion on dLl holds for 4 = 2, 3.' 


du 9V 


n n—l^l + l 

^-2^ _LgKI + |j|; 


(4.14) 


14=1 lil=i 

n+2-\i\-\j\-d 


r- 


1 + 


+0{e^+^), 


E 

p=i 


d 


^+^-^Qp){d^U{z))] M,,di—G{x,z) 




where the remainders are is dominated by 11^2(90) for some C inde¬ 

pendent of X € dLl. Here G{x,z) is the Dirichlet Green function, Mij, i,j E N'^, are 
the generalized polarization tensors, and Qp is the operator defined in ( where 
Vijk is defined, in this case, by 


(4.15) 


Pli = ^Mi,dlSn{di{^G){;z))\,,=, 

J • 


Theorem O can be proved in the exactly same manner as Theorem IT. We 
begin with Lemma |2.5|. Then the same arguments give us 


n n-|i|+l 

u{x) = U{x) - e‘^-^ Y, -e'*l+!^l(5*i4)(2)Mij5^G(x, z) + 0(e'^+^). 

14=1 lil=i 


From this we can get ( 1.141) as before. 


We conclude this paper by making a remark. The following formula is not exactly 
an asymptotic formula. However, since the formula is simple and has some potential 
applicability in solving the inverse conductivity problem, we make a record of it as 
a theorem. 
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Theorem 4.3 We have 


n n-jil + l 

(4.16) u{x) = H{x) + e'^-2 ^ ^ -e\^\+\^\d^H{z)MijCpT{x - z) + 0(e‘^+"), 

1*1=1 |i|=i 


where x G Hq o^nd term is dominated by C\\g\\i 2 (^QQ'^e‘^~^'^ for some C de¬ 

pending only on cq, and H is given in (2.8). 
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